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This paper describes algorithms for predicting end- to- end perform- 
ance measurements in private networks. Service characteristics such 
as end-to-end blocking and delay are calculated based on point-to- 
point traffic data and a network routing guide. These techniques 
have been incorporated in the Enhanced Network Administration 
System. The system is routinely used by AT&T Long Lines and 
operating company network administrators for private network de- 
sign and evaluation. 

I. INTRODUCTION 

An important procedure in network design is performance predic- 
tion. In a private network environment, a traffic engineer must rec- 
ommend a network design that satisfies the customer's performance 
requirements. A set of computer programs called the Enhanced Net- 
work Administration System (knads) has been designed for adminis- 
tering Enhanced Private Switched Communications Service (epscs) 
networks and Electronic Tandem Switching ifts) networks. The 
EN ads Network Service Evaluator (neteval) uses point-to-point 
traffic data and a network muting guide to provide detailed ser- 
vice characteristics of a network design. This paper summarizes the 
neteval algorithms. 

In addition to characterizing network service, neteval is used, in 
conjunction with other enads modules, 1 to ensure that the service 
characteristics chosen by the customer are achieved. The Network 
Synthesis (netsyn) module designs a network that is close to the 
service specified by the customer The service evaluator is used to 
ensure that the final network recommendation meets the customer's 
service requirements* 

Traditionally, performance estimates of a large network in an alter- 
nate routing environment have been confined to a segment of a 
network, such as the intermachine trunk portion. 2 neteval contains 
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Fig, 1 — Private network configuration 



decomposition and aggregation procedures that permit evaluation of 
complete end-to-end performance of a network as opposed to just a 
network segment. In addition, the enads ovaluator generalizes the 
Katz algorithm to include private network features such as queues, 
and controlled access to network facilities via Facility Restriction 
Levels (frls). This paper describes the decomposition and aggregation 
procedures used in the evaluator, as well as the additional tools 
required to analyze some of the network features germane to private 
networks. 

IK PRIVATE NETWORK CHARACTERISTICS 

A private network (see Fig. 1) interconnects customer locations f on- 
net service points) and other locations (off-net service points). An on- 
net service point is generally a PBX/Centrex or a key set connected 
directly to an access line group. Off-net points are served by off-net 
facilities (Foreign Exchange, WATS, or local off-net access lines). In 
this paper both access lines and off-net facilities are called end-links. 
One end of each end-link is associated with a service point and the 
other end is connected {or homed) to a switch. A network of interma- 
chine trunk groups connects the switches together, forming the inter- 
machine trunk-group portion of the network. Switches permit concen- 
tration of point-to-point traffic. 

Traffic generally originates at a service point, seizes an access line, 
and arrives at a switch, (It is possible in an ets network for a service 
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point to reside at the same location as the switch, so that an access 
line is not required.) neteval assumes that the switch examines only 
the final destination of the traffic and its frl and then associates an 
ordered route list of links at that switch. The route list is scanned until 
an idle circuit is found If no such circuit exists, the call is blocked (in 
the absence of queuing). If a circuit is seized, the traffic arrives at the 
next switch, and the routing procedure continues until the call is either 
completed or blocked. 

Other private network characteristics of interest for evaluation are 
queuing and frls. Queuing permits a call to wait at a specified trunk 
group for an idle circuit. The call is queued on a trunk group after a 
search through an ordered route list fails to find an idle circuit. The 
queued call may still be sent to reorder if the caller's wait exceeds a 
time-out threshold. Also, the caller may abandon the queue. The 
queue discipline is first come, first served* 

frls provide the capability to restrict or expand access to network* 
facility route lists. The frl can be based on the calling station and/or 
an authorization code. This provides the customer the opportunity to 
offer different grades of service by restrictive routing to different 
groups of users in the network. The frl and final destination of a call 
are used in deciding to which links in the network the call will have 



III. EARLIER SERVICE EVALUATION ALGORITHMS 

Statistical techniques are available for estimating service character- 
istics of networks. It is well known that if first-offered traffic is 
generated from a Poisson process and holding times are exponentially 
distributed, then network characteristics can be obtained by formulat- 
ing an appropriate Markov chain model and solving the resultant 
system of birth-dealh equations. This technique, however, requires an 
exorbitant amount of storage and computer time and is not considered 
feasible for large-scale networks* 

Another reasonably accurate statistical technique that requires 
knowledge of only the traffic mean and variance (or the variance-to* 
mean ratio, called peakedness) is the Katz algorithm/ which was 
originally designed for evaluating switch-to-switch blocking probabili- 
ties. This algorithm requires estimates of switch- to- switch traffic 
means and variances that are "effectively" offered to each trunk group 
in the network. 

Traffic that is carried on a link may be blocked on subsequent links. 
The holding time for such traffic on seized circuits is substantially less 
than that for completed traffic on the network. Effective* offered traffic 
reflects the shorter holding times of subsequently blocked traffic. The 
effective-offered load is used to compute the link-blocking probabilities 
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and traffic-overflow variance for each link. The algorithm is an itera- 
tive process that updates effective-offered loads and link parameters 
at the end of each iteration. 

After these initial link parameters have been computed, the switch- 
to-switch loads are distributed through the network based on the 
network -alternate routing plan and the link parameters. As this traffic 
is routed through the network, the effective means and variances of 
the traffic offered to each link are accumulated to update the link 
parameters for the subsequent iteration. The blocked traffic is also 
accumulated to provide probabilities for that iteration. After all the 
switch-to-switch loads have been distributed throughout the network, 
the link parameters are updated and the entire procedure is repeated 
until convergence in the swiich-to-switch loss probabilities is obtained. 

The Katz algorithm provides switch-to-switch blocking probabilities. 
Unfortunately, the number of calculations per iteration during the load 
assignment process is a function of the number of switches in the 
network. Traffic associated with each switch pair must be offered 
separately to each link in its routing path. Thus, an fr-switch network 
must distribute N{N - 1> switch-to-switch loads through the inter- 
machine trunk network. Since private networks generally consist of 
several hundred service points, it is not computationally feasible to 
define each service point as a switch and use the Katz algorithm. In 
addition, the Katz algorithm does not analyze queuing. 

IV. NETEVAL ALGORITHM 

neteval is an iterative algorithm that successively updates traffic 
loads and associated traffic characteristics until a convergence criterion 
is satisfied. Because of the size of the network, neteval decomposes 
the network into an end-link portion and an extended-trunk portion 
(see Section 4.4) to compute the traffic characteristics. After conver- 
gence the network performance components are aggregated to obtain 
end-to-end network performance estimates. 

Sections 4.1 and 4.2 contain the neteval assumptions and a brief 
outline of the algorithm. Subsequent sections explain in more detail 
the decomposition procedure, calculation of network parameters, and 
the aggregation techniques. 

4,1 Model assumptions 

The model assumptions for neteval are as follows: 
<i> Traffic means and variances provide a sufficient description of 
the loads. 

(it) The holding time on a link has four components: actual message 
time, ringing time, link set-up times, and subsequent queue delays. 
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Under this assumption blocked traffic can contribute positive loads lo 
the network, 

(iii) A blocked customer will redial with a specified retrial proba- 
bility and a variance that is a fraction of the overflow variance. 

1 1 ■■■ i No queue time-outs or abandonments occur. 

(v) The system is in statistical equilibrium. 

i vi ' Only the final destination and frl of a call are used at a switch 
tx> determine routing through the network. 

4.2 NETEVAL algorithm for service evaluation 

The basic algorithm procedure is given below. Some of the terms 
used are explained more clearly in later subsections. 

(i) Decompose the network based on routing into an end-link 
network portion and an extended trunk portion. 

Ui) Associate traffic parcels with each network portion. 
(iti) Initialize all parcel blockings* (Zero can be used if no other 
estimates are available.) 

(iv) Compute effective-offered parcels to the end-link network 
portion. 

id Calculate blocking probabilities for the end-link parcels. 
{vi) Compute total switch-to-switch and switch* to- final-destina- 
tion parcels offered to the extended trunk portion of the network. 
(vii) Calculate switch-to-switch and switch-to-final -destination 

blocking probabilities. 

(i;i*i) If end-link, switch-to-switch, and switch-to-final-destination 
blocking probabilities change significantly, return to step iv. Otherwise, 
go to step u 

(ix) Compute point-to-point characteristics* 

4.3 Decomposition 

As we mentioned in Section III, a typical private network is too 
large to model each service point as a switch and use the Katz 
algorithm. The network must be decomposed into segments and ana- 
lyzed separately. To define such a decomposition we must first define 
a traffic parcel. 

The aggregate of point-to-point loads with identical routing, when 
offered to a particular portion of the network, will be called a traffic 
parcel for that network segment. For example, if all point-to-point 
traffic originating on service points homed on switch / and destined 
for service points homed on switch J have identical route lists when 
offered to the inter-machine trunk-group portion of the network, then 
such traffic forms a switch / to switch J parcel in the inter-machine 
trunk-group subnetwork. Any technique using such a traffic aggrega- 
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tion assumes that the network characteristics of the traffic parcel 
sufficiently approximate those of the individual point-to-point loads. 

One method of decomposition is to separate the network into an 
inter-machine trunk-group portion and an end-link portion. Switch-to- 
switch parcels, as defined above, are associated with the inter-machine 
trunk-group network segment. All point-to-point traffic in the end-link 
segment, originating from a point to a switch and offered to the same 
route list of end-links, forms an originating end-link parcel. Similarly, 
all point-to-point traffic offered in the end-link network from a switch 
to a service point, using the same route list, is aggregated into a 
terminating end-link parcel. 

Since each parcel's point-to-point load components have identical 
routing in its associated network segment, blocking can be computed 
for each parcel using the Katz algorithm. For example, switch-to- 
switch blocking can be computed for switch-to-switch parcels in the 
inter-machine trunk-group network segment, and originating and ter- 
minating parcel blockings can be computed in the end-link segment. 

The above decomposition facilitates aggregation of parcel charac* 
teristics to obtain end-to-end characteristics. For example, if points i 
and j are homed on switches /and J, respectively, then i-to-y blocking 
(Mis 

bum 1 - (1- BU.HI - B U ){\ - B/i*rJ- (1) 

(See Appendix A for a summary of the notations used.) 

The three factors in the product represent the probability of call 
completion for the originating parcel i» the switch / to switch J parcel, 
and the terminating parcel y, respectively. An implicit assumption in 
such a decomposition is that all Mo*/ traffic must use end-links homed 
on switch J for call completion. Or equivalently, the last switch, called 
the terminating switch, that completed i-to-y traffic encounters must 
be switch J. Figure I, however, displays a typical private network in 
which terminating switches are not unique. 

It illustrates a two-level hierarchical trunk-structure, with on-net 
service points homed on lower level switches / and J t respectively. A 
bypass access line group is used to route traffic from switch K to point 
j. The two-level trunk hierarchy permits i-to-/ traffic to reach pointy 
using either switch J or K as a terminating switch* 

Since i-ta-j traffic is not required to use the home access line group 
serving j $ the above equation for 6 V is inaccurate. The last two terms 
in the product of the equation do not provide the switch / to point / 
blocking. To handle routing patterns with nonunique terminating 
switches, it is necessary to introduce a switch- to- final-destination 
parcel. Such a parcel represents all point-to-point loads from service 
points homed on a switch that has identical routing from the switch to 
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the destination service point. In Fig- 1, a switch / to point ; parcel is 
required to compute Uta-j blocking. If B° is the blocking probability 
for such a parcel on the network, then b t , can be accurately estimated 
by 

*6- l-U-Bj^Hl-B*). <2> 

Thus, the neteval decomposition must classify traffic into four 
parcels for adequate performance estimation; originating end-link par- 
cels, terminating end-link parcels, switch-to-switch parcels, and switch- 
to- final-destination parcels. Switch* to- final-destination parcels are 
used only for service points that do not have unique terminating 
switches. To compute characteristics for switch-to-switch and switch- 
to-final-destination parcels, a segment of the network must be defined 
that contains all links used by these parcels. This subnetwork consists 
of all inter-machine trunk groups and those end-links contained in the 
route lists of the switch- to-final-destination parcels. Such a collection 
of links is called an extended trunk network in this paper. With the 
extended trunk network and its associated switch-to-switch and 
switch-to-final-destination parcel means and variances, the appropri- 
ate parcel characteristics can be computed. Similarly, the end-link 
subnetwork and end-link parcels permit the calculation of end-link 
parcel characteristics. Section 4,6 describes how parcel characteristics 
can be aggregated to form the desired end-to-end characteristics. 

4.4 Link analysis 

Both the end-link and switch-to-switch and switch-to-final-destina- 
tion analyses are based on offering effective- offered loads to a single 
trunk group and computing associated link characteristics that are a 
function of overflow parcel means and variances. Whenever alternate 
routing is involved in a network segment, the Katz algorithm is used 
on the appropriate section. (Even in a queuing environment this 
procedure is valid except that a new algorithm to compute link param- 
eters is used,) frl analysis is included in neteval by stratifying parcel 
loads by FRL and associating different route lists with each frl 

grouping. 

Effective-offered loads are computed for each link, from which link 
parameters such as blocking and overflow variance are estimated. The 
effective-offered loads associated with a link must take into account 
interactions with the rest of the network. Traffic blocking both prior 
and subsequent to a parcel being offered to a link reduces the effective- 
offered traffic load to the link. However, network setup times, ring 
times, queue delays, and retrial attempts increase effective- offered link 
loads. Even though the end-link and extended trunk network segments 
are analyzed separately in neteval. the effective-offered load equa- 
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Fig. 2 — Link configuration. 

tions are designed to reflect the above network characteristics and 
thus account for the interactions between the network segments. These 
interactions are updated successively in the form of updated parcels in 
steps iv and vi as summarized in Section 4-2. 

The following subsections describe the computation of link param- 
eters both with and without queues. 

4.4.1 Unqueuedcase 

Figure 2 displays a general link configuration with link-offered 
parcels (a,, ft), i = 1 to m, a, is the mean of the *th parcel and v t is its 
variance. [If queues are present, the queue parcels are (a tff , t-v).] The 
equivalent random method can be used to compute overflow means 
and variances for aggregate traffic (a, u). To apportion the aggregate 
overflow moments to different parcels we use the Katz parcel -splitting 
approach. If 6, is the e'th parcel blocking, then 

where, 

b - aggregate link blocking, 

z Q - peakedness of aggregate traffic, 

z, - Vifa tt and 

k « modification factor. 

Several heuristic formulas can be used for A, Katz 2 expresses A as a 
function of z and h . We have selected the heuristic formula from the 
Defense Communications Engineering Center 

k - «***, (4) 

where c - 2$*QzZ zm t t - a0528*; 4l€3 , s m 5A56z; 2Am t and n is the 
number of trunks. 

Once 6, is known, the parcel overflow and carried mean estimates 
are a,bi and a, (1 - 6, >, respectively. The overflow and carried variance 
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are approximated by £*<*,&■ and u, — .".■■'■. where zt is the overall 
overflow peakedness. The overflow variance estimate results in equal 
peakedness to all parcels. 

4.4.2 Queued case 

To derive the overflow mean and variance for a parcel offered to a 
link in the presence of queues, another approach is needed. We 
approximate the behavior of a link containing queues with a two- 
dimensional Markov chain for each peaked parcel and a one-dimen- 
sional Markov chain for each smooth parcel The parcel is modeled 
based on its mean and variance. If the parcel is peaked (peakedness 
greater than one), an Interrupted Poisson Process <ipp) is used. Smooth 
parcels (peakedness less than one), which can MOir when traffic 
carried on previous groups is offered to a trunk group, are modeled as 
a Poisson process with parameters adjusted to better reflect the small 

peakedness. 

4.4.2.1 Interrupted Poisson processes. An ipp. also called switched 

Poisson, is a Poisson process with rate K which alternately for some 
period of time shuts off all arrivals from the Poisson process, and lets 
the arrivals go through for another period of time. Both time periods 
are exponentially distributed, independent of each other and any 
previous time periods, ipp is used to model traffic with peakedness of 
one or greater. We say that the switch is ON when the arrivals go 
through, and otherwise the switch is OFF- The IPP process is shown in 
Fig. 3. A Poisson process is an ipp with the switch constantly in the 
ON position. Given (a,, v,) t an ipp with an expected on time y~\ an 
expected off time w" 1 , and a Poisson rate A (when the ipp is in the ON 
state) can be obtained from the relationships 

and (5) 



L_' 



+ A 



it, - a, + of - a,\ — , , (6) 
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Fig- 4 — Markov chain for peaked parcel inpul. 

where A is arbitrarily set to the larger of a, and v t . Vt — a, + a? is the 
second factorial moment of the m\ 

This completely specifies the ipp point process. For further proper- 
ties of the ipp, see Kuczura. 3 

4.4.2.2 Peaked parcel analysis. The ipp ON and OFF states form one 
dimension of the Markov chain used in the queued case. The other 
dimension consists of three superstates: 
(i) Not-all-trunks-busy 
(it) All- trunks- busy, but no-calls-in-queue 

(m) Some-calls- waiting-in-queue. 
Note that in case (;/n. necessarily all trunks are busy, and that the 
system moves from state 1 only to state 2. from state 2 either to state 
1 or to state 3, and from state 3 only to state 2. No transition between 
states I and 3 is possible. 

Thus, the Markov chain space is (X t V). X - 0, I, Y - 1, 2 t 3, where 
X is the ipp status of a parcel offered to a link and Y is one of the three 
superstates discussed above. Figure 4 displays the flows among the 
states. 

Modeling the complex behavior of the trunk group with queues by 
compressing all states into superstates and assuming constant transi- 
tion rates between them is a crucial assumption. The simplified state 
space is not Markovian any more, and the rates between the states are 
generally nonconstant. However, for first-moment calculations, such 
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as the mean of the overflow stream, the replacement of the noncon- 
stant transition rate by the average rate gives exact results. The 
calculation of the rates is given below. Some of the rates are themselves 
approximations for the average rate. 

Since traffic is offered to the queues only when all trunks are busy 
in the group, we will consider the queue-offered traffic "conditional on 
all trunks busy." and use the conditional mean and variance of these 
streams when computing link characteristics. If (a v , v<,) is the uncon- 
ditional queue mean and variance, and p is the probability that all 
trunks are busy, then the mean and variance conditional on all-trunks- 
busy is given by 

These equations are based on the assumption that if the random 
variable X is the number of busy servers in an infinite trunk group 
receiving traffic from the queue stream, then EX* = pEX' Ct i - 1, 2. 

4.4.2.3 Transition rates. To compute transition rates w? N and &>i , 
we need to define u- as the average frequency of transition from not- 
all- trunks- busy states to all- trunks-busy states. Since this rate does 
not depend on the presence of queues, we may calculate wi on a trunk 
group without queue*. I-et h T be the time congestion of a trunk group 
of size n, offered traffic with mean a and variance v. Further, let T\ 
and T-2 be the expected sojourn times of the trunk-group system 
(without queues) in states not -all- trunks-busy and all- trunks- busy, 

respectively- In that case 

Since Ta • n~ l (assuming unit average holding time) for known ftr, 
one may find 7\.and hence u?i = \/T u immediately. The time conges- 
tion b : is equal to the blocking, b t , that would be experienced by a 
Poisson parcel offered to the trunk group. This is given by the Katz 
parcel splitting formula 

6, -Ml + *<*-**>], <10) 

where z, is equal to 1. 

Now that ui is known we can solve for <j? N and w? FF . We require 

that in the long run the number of transitions from states [(1, 1), 

(0, 1)] to [(1, 2), {0, 2)] in Fig. 4 is u?i per unit of time. Since movement 

among these states does not involve transition rales in or out of states 

[(1, 3), (0, 3)1 it is more convenient to consider a smaller Markov chain 
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(X\ y*>,X* =0,1, y* = 1, 2, which is a subset ofthe chain shown in 



OFF „j ON 



Fig. 4. We can then express the relationship between &}, wV , and <j\ 



« , P 01 „OFF , P" „ON / in 

Poi + pu P01 + Pll 

where p v is the steady-state probability distribution for the Markov 
chain (X\ Y m ). 

We still require an additional constraint to make w? FF and w? N 
unique. A reasonable additional assumption is that in the absence of 
queuing, the parcel blocking is the same as that obtained from the 
parcel splitting formula. Thus, we require 

4. *!—. (12) 

Pll +P12 

where b, is obtained from the Katz parcel splitting formula* The p<, are 
themselves expressed in terms of «r* and wV N . However, it can be 
shown that the equations above lead to an equation no worse than 
quadratic. 

To compute the transition rate y it assume the system is in an all- 
trunks-busy state, with no-calls-in-queue. yi is the rate with which 
trunks become available. Since we express time in multiples of an 
average holding time, yi = n (the number of trunks). Further, under 
the assumption that the holding times are all independent, exponen- 
tially distributed variables, the transition rate .■ . is a constant, inde- 
pendent of the time spent in the all-trunks-busy state, with no-calls- 
in-queue. 

To compute uz assume the system is in the state with all-trunks- 
busy and no-calls-in-queue. Then the frequency with which the state 
some-calls-waiting-in-queue is entered is the input rate into the queues. 
Since the queue-input rates a . and n... are given conditionally, it is 
clear that 

W* ■ acq t + Oc«|i (13) 

where u>2 is a constant, rather than the average rate when the queue 
input streams are both Poisson. If the input streams are peaked, they 
may be represented as interrupted Poisson processes, and a more 
precise, but more elaborate calculation for <*>i may be carried out, 
taking account of the variances of the queue-input streams. This 
precision was not considered necessary* 

Note that in all of the above, it is possible that one or both queues 
are non-existent. The mean Oo? should be set equal to for non-existent 
queues, 

yz, the rate from the state call-in-queue to the state with all-trunks- 
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busy and no-calls-in-queue, is the inverse of the mean length of time 
the system spends in the state with some-calls-in-queue. What we 
therefore need to calculate is the average duration from the time the 
first call enters one of the queues until the first subsequent time both 
queues are empty. We denote this mean duration "by Q and call it "the 
queues busy period/' One then has 

y, - Q- 1 . (14) 

The calculation of Q is appreciably simpler for links with one queue 
than for links with two queues* In either case, one models the input 
stream(s) into the queue(s) as ipp(s), and then finds the simultaneous 
ergodic distribution off the number of waiting calls and the ipp input 
state. This ergodic distribution is a conditional distribution, given that 
all trunks are busy- The input streams are taken to be conditional on 

this event* 

Lastly, cj,> and y* are the on and off rates for an ipp parcel offered to 
a trunk group- Thus, eqs. (5) and (6) can be used to compute w and 

4.4,2.4 Computation of link characteristics. The overflow mean, 
overflow variance, and carried variance of a parcel <a, , v t ) are computed 
by adding another dimension to the state space, which is the number 
of calls in an infinite trunk group that receives parcel overflow traffic 
when the system is in an all- trunks- busy state. Moment-generating 
functions can be uaed to compute the deaired quantities. (Details are 
provided in Appendix BO The overall link blocking is then the ratio of 
the sum of the parcel overflow means to the sum of the parcel offered 
loads. 

Bach of the parcels offered to the queues gives rise to two resultant 
parcels, namely the overflow parcel and the carried parcel. The queue 
overflow parcel consists of those calls finding all queue slots occupied, 
and the carried parcel contains all other calls, since we assume that 
there are no abandonments or time-outa. 

As a result of assumed independence of the parcels offered to the 
group, the parcels offered to the two queues are independent during 
the time intervals that all trunks are busy. This means that waiting 
times and other variables associated with a two-queue trunk group 
may be estimated through simple Markov chain calculations. 

Queue blocking and delay are computed conditional on all-trunks- 
busy, A Markov chain (/i, ju *i. h), where the y^s are the number in 
each queue and the i*u represent the on/off conditions of IPP queue 
streams, is used to compute these quantities. (Appendix C explains the 
algorithms used to compute the steady-state probabilities of the chain,) 
For example, if e /ithJlJi is the ergodic distribution of the chain and 
(yi* w>) represents the on and off rates of the Q* arrival stream, then 
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Fig. 5 — Markov chain for smooth parcel inpul. 

the Qi blocking probability is 

&,, -fl+Iij Jt^io + e^ii). (15) 

Since a customer can arrive at Qi only when the ipp is ON, eq. (15) 
is expressed only if the queue arrival stream turned on* 

4A.2.S Smooth traffic parcel analysis. A smooth traffic parcel is 
offered to a link when v t is less than a, . A somewha t simplistic approach 
is used to compute the overflow stream characteristic for a smooth 
parcel A one-dimensional Markov chain is formed, assuming the 
parcel is Poisson (see Fig, 5|. All rates, except for wf M (used in place of 



Ul 



ON „j OFF 



and «i ), are the same as in the case of the peaked traffic analysis. 
o>? M is computed in the same manner as oj\ , but with zt set to the parcel 
peakedness. 

Overflow quantities can be calculated using the same techniques 
applied to peaked parcels. The corresponding variance calculation, 
however, leads to too large a variance of the overflow, because the 
input into the trunk-group characterization switch is taken to be a 
Poisson stream with mean and variance (a lt a), instead of the given 
smooth stream with mean and variance (a, t w). To correct for this 
effect, we reduce the calculated variance by multiplying it by the 
peakedness (less than one), of the offered parcel. The estimated 
overflow variance is z r v?* t where vf* is the overflow variance based on 
Poisson input* The mean overflow is given by 

m cof M (^ + n) 
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4.5 Convergence criteria 

The convergence criteria for step viii of the neteval algorithm 
(Section 4.2) assumes convergence has occurred if the absolute differ- 
ence of each end-link, switch-to-switch, and switch- to- final-destination 
blocking probability in successive iterations is within a specified limit. 

4.6 Point-to-point characteristics 

Once the network has been decomposed into the extended trunk 
network and the end-link network, and the components have been 
analyzed, aggregation techniques are required to compute point-to- 
point characteristics. This is done by associating two or three traffic 
parcels, depending on the routing patterns, with each point-to-point 
pair. If *-to-y traffic is part of a switch- to-final -destination parcel when 
offered to the extended trunk network, then only two parcels are 
associated with the point pair: the originating parcel at point i offered 
to the end-link network segment, and the switch-to-final-destination 
parcel in which it is contained. If the i-to-y traffic is part of a switch- 
to-switch parcel in the extended network, then three traffic parcels are 
associated with the point pair: the originating parcel, the switch-to- 
switch parcel, and the terminating parcel in which i-to-y traffic is 
contained. Point-to-point characteristics are functions of the corre- 
sponding parcel characteristics. 

If i*tO-7 traffic is part of a switch-to-switch parcel when traversing 
the extended trunk network, the i-to-y blocking by is 

6„ = 1 - (1 - BlomJU - B u )i\ - B/J. (17) 

If i-tO-y traffic, however, is contained in a switch-to-final-destination 
parcel, 

6t,-l-(l-Bioct>U-B*), (18) 

The i-to-y expected queue delay, W tJt is a weighted average of the 
expected queue delay on trunks and end-links. When i-to-j traffic is 
contained in a switch-to-switch parcel, 

W« - (1 - BlmjW + (1 - BU)d - B u )Wf. (19) 

The trunk delay, W u % and terminating link delay, Wf % are weighted 
by the probability that e-to-y traffic will be offered to a link on which 
it is queue-eligible. When j-to-/ traffic is contained in a switch-to-final- 
destination parcel* 

where W h is the expected queue delay for traffic originating at switch 
/ destined for point j. 1 - B' l<>n is the probability that the traffic seizes 
an end-link and arrives at switch J. 
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Fig. fc— End-to-end blocking distribution— Network A. 

The i-to-y queue delay probabilities, P/) VI are computed analogous 
to Wiy. When i-to-y traffic is contained in a switch -to-switch parcel 
when traversing the extended trunk network. 

PD V - (1 - BU n )PD u + (1 - BU)d - B"-PD tJ )PDf t (21) 

where (1 - B' WK )(1 - B /,; - PZ) w )PO/, is the probability of traversing 
the trunk portion of the network without queue delay but incurring 
queue delay on the terminating end-link. 

If t-to-y traffic is contained in a switch -to* final-destination parcel. 

PD V - (1 - Bf WH )PD°. (22) 

Conditional expected queue delays are given by Wg/PD^t for PD V > 
0. Note that by creating an artificial access-line group with no blocking 
or delay, the above equations remain valid for ets networks with 
traffic originating at a tandem switch. 

V. APPLICATIONS 

epscs and ets networks generate call-detail records- These records 
contain information on individual calls in the network from which it is 
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possible to estimate the end-to-end blocking distribution in existing 
networks. This information can be compared with an en ads neteval 
analysis of the same network. Figures 6 and 7 display both call detail 
and netkv a L- derived blocking distributions for two existing networks. 
The means and variances in Fig. 6 are very close. The variances in 
Fig* 7 are not as similar as their corresponding means. This is primarily 
because the midpoint of the last cell in the call-deuil-derived distri- 
bution is 43 percent, compared with 30 percent in the neteval 
distribution. Note, however, that the overall shape of the distributions 
within Figs* 6 and 7 is similar* This suggests that the modeling 
techniques in our evaluation algorithm provide a reasonable prediction 
of end-to-end performance. 



VI. SUMMARY 

Algorithms for computing end-to-end and link characteristics in 
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private networks have been described in this paper The algorithms 
generalize the Kata procedure to provide end-to-end characteristics for 
large-scale networks and also model private network features such as 
queuing and frls. These techniques have been incorporated into the 
enads Service Evaluator. neteval performance predictions have 
compared well with actual measurements for several networks. The 
enads neteval module is now routinely used with the companion 
netsyn module by AT&T Long Lines and Operating Company private 
network administrators for network design and evaluation. 
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APPENDIX A 

Notations 

Much of the notation used in Section IV is summarized below for 
the reader's convenience. 
Parcel Characteristics 
Blocking 

b v * Point i to point / blocking. 
B'oir ■ Originating end-link blocking for service point i homed 
on switch /. 
Bf itrm ■ Terminating end-link blocking for service point / homed 
on switch J. 
B u ■ Switch / to switch J blocking. 
B iJ — Switch / to final destination point / blocking. 
Delays 

Wy ■ Point i to point j queue delay. 

Wf = Terminating end-link queue delay for service point / 

homed on .switch J. 
W u = Switch / to switch J queue delay. 
W 1 ' = Switch / to final destination point j queue delay. 
PDjj = Point j to point j delay probability. 
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PDJ = Terminating end-link delay probability for service point 

/ homed on switch J. 
PD rj = Switch / to switch J delay probability. 
PD h = Switch / to final destination point j delay probability. 
Link Parameters 

(a,, v t ) - Mean and variance of eth parcel offered to trunk group. 
z, ■■ Peakedness of jth parcel offered to trunk group, 
(a 9it v qi ) — Mean and variance of load offered to queue. 
(a fV(l i\ Vi ) — Mean and variance of load offered to queue conditioned 
on all trunks busy. 
bi = Blocking of ith parcel offered to trunk group. 
(^ u) = Aggregate trunk-group load mean and variance. 
z* = Trunk-group peakedness of offered load. 
zt m Trunk-group overflow peakedness. 
b - Trunk-group blocking. 
Markov Transition Rates in Trunk-Group Queuing Model 
y„ = Interrupted Poisson Process off rate- 
id,, = Interrupted Poisson Process on rate. 
u? N - Not-all-trunks-busy to all- trunks- busy, no-calls-in-queue 
flow rate t given Interrupted Poisson Process (ipp) on. 
c*)2- All -trunks-busy, no-calls-in-queue to calls-waiting-in- 
queue flow rate. 
Jf n = Not-all-trunks-busy to all-tninks-busy, no-calls-in-queue 
flow rate, given ipp off. 
y* ■ Calls * waiting -in -queue to all -trunks -busy, no-calls-in- 
queue flow rate, 
y, - All- trunks-busy, no-calls-in-queue to not -all -trunks-busy 
flow rate. 

APPENDIX B 

Overflow Mean and Variance of Peaked Parcel 

This appendix gives the details of the calculation of the overflow 
mean and variance of a parcel, modeled as an IPP, ottered to a trunk 
group, which in turn is modeled by three states and the four transition 
rates among them. Let the variables \ yo, yi, ya< <**• W ? N ' W ? FF anc * U2 
be as in Section 4.4.2. 

The system to be treated, including the infinite trunk group, is a 
continuous- time Markov chain, with state space 

N« (<A, i;»:* = 0, 1, .- ;i-0, 1;>= 1,2,3), 

where k is the number of calls in the infinite trunk group, i is the status 
(ON - 1, OFF = 0) of the ipp, and / is the state of the trunk group, 
(not-all- trunks-busy, all-trunks- busy, but no-calls-in-queue, some- 
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calls* waiting- in queue). Let > denote the status of the jit and S? the 
status of the trunk group. Designating by p* tJ the steady-state proba- 
bility that the system is in state (A, i, /), we define the generating 
functions 

■»(■)- 2 Pm»*t for 055S1, i = 0,i, y-1,2,3. 

The mean number of calls in the infinite trunk group is 

ui k u as as 

where r stands for the six element row vector of ones and 



;is) 






.Similarly* the second factorial moment is 

EX(X - 1) = 1 1 k(k - Dpm, - I -^ v v {l) - 1'-£ *w> 

from which one easily finds the variance. Writing out the steady-state 
equations for pk% } and summing appropriately, one finds that the ir v (s) 
satisfy a system of equations, which we will define in vector form. We 
define 












tff rK (uu + ri + «*) >J 

tfe -((*> + ?*) U 

<* o o -(y* + «*n yi o 

«* w?" -(>* + y t + un) yi 

«o *» -!)*+?*) 



and the diagonal matrix B by 

B-diag(0,0 ( 0,0,A,A), 

where the elements of the diagonal are given inside the parentheses. 
The vector n(s) satisfies the equation 



is - 1) — w(s> - Avis) + (s - l)Bv{sh 
ds 



(23) 



Setting s= 1 in (23), one obtains Am(1 ) - 0, from which the marginal 
probabilities ir v (l) = P{Si » e, S? — J) are found. Taking derivatives 



2162 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1982 



wj.t & in (23) and setting & - 1 gives 



CM 

Since l'/l =0, one gets 

i'4-bQ) - n/- aj^»(d - L'Bffd), (24) 

and hence the mean a * of the overflow may be found by solving for 
Ait(\) - 0, with the &ide condition Vw(l) - 1. Differentiating eq. (23) 
twice, one gets 

+ 2B — *(«) + (a - 1)B — , *<*>- 

ds as 



d' 



Setting 8=1 and solving for— ^ tt(1) gives 

* 7r(l(=2(2/-A> l B — »q). 
tfs 2 ds 

Again, since l'(2/ - A) - 2 1'. one finds 

l'iir(l)- I'B-^wUh 
ds* ds 

or, substituting (24), 

ds* 
Thus, by inverting the 6x6 matrix I — A% the variance of the overflow 

may be found* 
The carried traffic mean may be obtained directly from the overflow 

mean via 

a r - a- a*> 

where a t is the carried mean, a the offered traffic mean, and a' the 
overflow mean. It can be shown that a? may alternatively be found via 
the procedure sketched above, if we change the B matrix to 

B-diag(0,0,0,A,0,0). 

Thin change represents the fact that carried traffic leaves when the 
ii'p switch is ON, and the trunk group characterization is in the not- 
all- trunks- busy state. 

With this change in B. the variance of the carried traffic may be 
calculated in exactly the same way as the overflow variance. 
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APPENDIX C 

Queue Analysis 

This appendix gives some details of the queue analysis. Assume a 
trunk group with n trunks, q% queue slots on one side, qi on the other, 
and input streams into the queues given, conditionally on all-trunks- 
busy, by their ipp parameters (Af, <*■ Yih i ■ 1» 2. The system is a 
continuous time Markov chain. This appendix gives a description and 
a fixed ordering of the states of the system in Section C.l and 
algorithms for the steady-state probabilities in Section C.2. 

C. 1 Ordering of the states 

The elements of the different vectors and matrices to be used in this 
appendix are characterized by a four-dimensionaJ index, (*, j % ol, o2), 
where each possible index is a state of the system* Index i is the 
number of calls waiting in the first queue, j the number of calls waiting 
in the second queue, and ol and o2 are status bits, indicating whether 
the ipp representing the input stream is currently in the OFF state 
(o = 0) or in the ON state (o - 1). Status bit ol refers to the first 
queue. o2 to the second. Index j runs between and q Xt J between 
and q 2 . All bounds are inclusive, 

A fixed order of enumeration is adhered to. There is a macro 
ordering, for the (i, j) part of the state designations, and within that 
ordering, the status bits have a fixed order. 

The gross order is given by 

(2,0), ..- (q,,0), 

iqu 1), 



(0, 0), 


(1,0), 


(0, 1), 


(1. 1), 


(0, 2), 


• 



<0 (9a ), (1,?*), 



(quqil 



Within these index elements there are four states indicating the 
status bits, which are ordered like 00, 01, 10, 11. Thus, the overall order 
is indicated by 0000, 0001,0010, 0011, 1000, 1001, 1010, 1011, 2000, etc 

C.2 Steady-State Equations 

The steady-state equations form a set of 4 (q t + I) (q^ + I) equations 
in as many unknowns, with the normalizing side conditions that the 
probabilities add to 1. In the ordering given above, the structure of 
these equations is that of a tri-diagonal blocked matrix equation, as 
follows. Let € f be the vector of ergodic probabilities associated with 
the states having j calls in Qi* j m 0, 1, ■ • • f qz* i.e., 

*y ~ (ffp/no, eo/oi. £o>io» ftyn ( eiyOH • • * * e«,>ii)- 
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Then the e, satisfies 

cw4, + e } C - 0, (25a) 

e>-ifl + e,Ae + ej*iC = 0, J« 1, ... 9 g*~ ] (25b) 

e Vj ~iB + e*.-4/-°- < 25c > 

Writing out these equations in matrix form will show the tri^diagonal 
block structure. 

The matrices A* t A/, A fit B, and C arc all square of order 4(</i + 1); 
in addition, B and C are both diagonal matrices. Below we discuss each 
of these matrices, and exhibit their structures. 

(i) C is associated with a call arriving into Q 2 : the system moves 
from state (i,j t OuOi) to state (*»/+ 1, oi, o 2 ), with transition rate fa 
if and only if 02 - L Therefore, C - diag(0, A 2 , 0, A 2 , • ■ ■ , 0, A 3 ). 

07) B signifies the transitions caused by a call waiting in Q 2 being 
served. The transition is from UJ t Qt*<h) to UJ- 1, o tt o 2 ), with rate 
n if* = 0, and with rate n 2 (where n 2 is the effective number of circuits 
available to calls from Q 2 when Qi and Q 2 are both non-empty) if i * 
0. Specifically, B - diag(n, n, n t n, Ra, rt2, * * • , na). 

(m) A-. -4/, and A* are all associated with transitions for which the 
number of calls in Qz is constant, this constant being equal to for A, 
(empty), q% for A f (full), and general for 4* , These matrices have a tri* 
diagonal block structure themselves, with blocks of size 4X4, with the 
super- and sub-diagonal blocks a&sociated with calls entering Qi and 
leaving Qi, respectively, and with the diagonal blocks modeling the 
transitions between the input status slates (oi, o 3 ) ■ (0, 0), (0, 1), (1, 
0), (1, 1). -4,, A f and A* have slight differences owing to the fact that 

(a) The diagonal element for any state contains the negative sum of 
the transition rates out of that state. 

(6) The rate at which calls leave Qi depends on whether Q* is 

empty- 

Next we show that we may take advantage of the block structure of 
the 4(<7i + IHq* + 1> set of steady-state equations, and reduce it to a 
set of size Mqi + 1) to be solved via brute-force matrix inversion 
methods. Consider equations (25). Since B is a nonsingular diagonal 
matrix, we may write 

e^-i=-e 92 tCB ' + A f B"A B-% 

In general, if e, = e^Dj fory - k + 1, * + 2, - - - . 92, one obtains D* 
from equation (25b) 

D* = -Dk+iA M b- 1 - D*+2CB ] (* - 0, 1, • • • f fli - 1), 

with D Ut = /, the unit matrix, and J3 V? *» " °- 
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Assume then, that we have expressed e in terms of e^ . Now we use 
the top equation to solve for e . and get 

«h(Zhl + i>iC) - 0. {26} 

The requirement that all probabilities sum to 1 can likewise be 
expressed in terms of e Vi only: 

Xe,l -Za^Dyl = L 

Thus, we sum the rowsums of the D /t obtaining 

X-2D/L 

We now replace one column of DoA e + D\C by the vector x, and the 
corresponding entry in the vector in (26) by 1, and solve this modified 
set of equations. We obtain «■.,., and from it may find all ergodic 
probabilities by postmultiplying e^ by the D/s* 

If the number of queue slots for Ql is larger than for Q2 (q\ > 92), 
we may choose to relabel Ql and Q2, and reduce the problem size by 
performing the recursive procedure described above over the larger of 
</i and </2- Thus, we may reduce the size of the matrix to be inverted 
from 4(<7i + l)(q 2 + 1) to 4[1 + min(?i, n)l 



2166 THE BELL SYSTEM TECHNICAL JOURNAL. NOVEMBER 1982 



